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Over a large fraction of pliase space a combination of an operator product 
and heavy quark expansions effectively turn the decay B -^ D^*'^e^e~ into 



PsJ ' a "short distance" process, i.e., one in which the weak and electromagnetic 
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^^ ' interactions occur through single local operators. These processes have an 

\^ • underlying W-exchange quark diagram topology and are therefore Cabibbo 
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On ' allowed but suppressed by combinatoric factors and short distance QCD cor- 

,^ ' rections. Our technique allows a clearer exploration of these effects. For 

o ■ the decay Bd^g — *■ J/'4'iVc)e~^e~ one must use a non-relativistic (NRQCD) 

T^ • expansion, in addition to an operator product expansion and a heavy quark 

•^ ' effective theory expansion. We estimate the decay rates for Bd^s -^ J/tpe'^e~ , 
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I. INTRODUCTION 

In a recent paper [|l[ we considered the collection of decays B^ -^ D^*^ e+e~. The decay 
rate for these is proportional to |Ktfep. We found that over a large kinematic domain one can 
reliably estimate the rate (in terms of |KibP)- The process is first order weak and first order 
electromagnetic, and, therefore, the amplitude involves long distance physics. The central 
observation of [|I| is that over a large kinematic domain the interaction is local on the scale 
of strong dynamics. The amplitude can, therefore, be approximated by the matrix elements 
of local operators, which can be estimated in a variety of ways and should eventually be 
determined in numerical simulations of QCD on the lattice. The branching fraction for 
B^ -^ D*^e^e~ , restricted to invariant mass of the e+e^ pair in excess of 1.0 GeV, was 
estimated to be 1.9 x 10~^. This is too small to be measured in e^e^ B-factories, but could 
be observable at high luminosity high energy hadronic colliders. 

In this paper we consider the decays Bgd — ^ J/4'e'^e^ , Bgd — >• 77^6+6^, Bs^d -^ D*^e^e~ 
and Bs^d -^ D^e^e^. These proceed via W-exchange topologies, as shown in Fig. |l[ In 
addition, 5^ ,, -^ J/ipe^e~ and Bd^s -^ VcG~^^~ have small contributions from penguins, 
which we neglect. The goal of the paper is to show how the methods introduced in paper 
|jl|] can be applied to the processes considered here. The kinematics of B^^s ~^ D^*-'°e^e~ is 
similar to that of B^ -^ Dj'l e+e~ so one expects the methods to apply readily. In fact, the 
only dynamical difference is that in Bd^s -^ D^*^^e^e~ the heavy b quark decays to a heavy 
c-quark, whereas in 5+ ^ D^ j e+e~ it is a heavy 6-anti-quark that decays into a heavy 
c-quark. The case Bd^s -^ J/'ipiVc)^'^^' is clearly different: both quark and anti-quark in 
the final state are heavy and they are moving together in a bound charmonium state. As 
we will see the expansion that arises naturally corresponds to NRQCD, the non-relativistic 
limit of heavy quarks bound by QCD into quarkonia. 

m ( ^^ I 

The processes under consideration here have advantages compared to B^ ^ D).^ e^e^. 
These processes are not suppressed by the small CKM element |\4fep. One might hope that 
the decay rate is, therefore, substantially higher. However, the enhancement of the rate due 
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FIG. 1. W-exchange quark topology diagram underlying the transition B^s — > D^*'^e~^e~. 
Emission of a e~^e~ pair from any line is understood. 

to bigger CKM elements is partially cancelled by small Wilson coefficients. Therefore, all 
these processes have small branching fractions. While none are observable at B-factories, 
some are observable at future hadronic collider experiments like LHC-B and BTeV. 

These processes are first order weak and first order electromagnetic, and, therefore, the 
amplitude involves long distance physics. We will show that over a large kinematic domain 
the interaction is approximated by a set of matrix elements of local operators. All these 
matrix elements should eventually be determined by lattice calculations. For the processes 
considered in this paper, the number of independent matrix elements is reduced by the use 
of rotational, heavy quark spin and chiral symmetries. 

This paper is organized as follows. In Sec. |I| we review the methods of Ref. [|l| that lead 
to an expansion in local operators. The review is done in terms of the graphs relevant to 
B -^ D^e^e^ , which is one of the processes of interest here. In Sec. |T| we present a novel 
analysis that shows that the matrix elements of the operators in the expansion are all related 
by a combination of heavy-spin, rotational and chiral symmetries. We then proceed to find 
the short distance QCD corrections to our operator expansion in Sec. |V]. In Sec. |V| we give 
expressions for the differential decay rates in terms of matrix elements of local operators. 
These should be considered our main results. To get some numerical estimates of the decay 
rates we crudely approximate the local matrix elements. The material in Sees. 0-0 deals 
with the decays Bg -^ D*^e'^e~ and Bq -^ D'^e~^e~, and we repeat the steps applied to 




u 



FIG. 2. Feynman diagram representing a contribution to the Green function. The filled square 
represents the four quark operator O and the cross represents the electromagnetic current jj^^, cf. 
Eq. (^), which here couples to the c-quark. 



the processes Bq -^ rjcC^e and Bg -^ J/ipe^e in Sec. Vl. Our results are summarized in 



Sec. VII. 



II. OPERATOR EXPANSION 

In this section we review the method introduced in [Q. However, we will present the 
method as applied to the process Bd — > D^*'>^e^e'~ . Therefore we will at once review the 
method and perform the necessary calculation for one of the cases of interest. 

The effective Hamiltonian for the weak transition in B^ — > D^*^^e^e~ ^ is 

AC 
Ke = -^ V^dV:, {c{fi/Mw)0 + cs{fi/Mw)Os) , (1) 

where 

= dYP-b c^uP-U (2) 

and 

Og = dYP-T% c^^P_T^u, (3) 

P± = (1 ± 75)/2 and T" are the generators of color gauge symmetry. This is a useful basis 
of operators for our purposes since the hadronic matrix element of the "octet" operator 
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FIG. 3. Same as Fig. g but with the electromagnetic current coupling to the 6-quark. 

Os is suppressed. The dependence on the renormalization point fi of the short distance 
coefficients c and Cg cancels the /i-dependence of operators, so matrix elements of the effective 
Hamiltonian are /i- independent. 

The amplitude for Bd -^ D^*''^ e'^ e^ , to leading order in weak and electromagnetic inter- 
actions and to all orders in the strong interactions involves the following non-local matrix 
element: 

{D*^\jd'x e^^- T{3Ux)0{mB^). (4) 

Here q denotes the momentum of the e^e~ pair, j^^ is the electromagnetic current operator 
and the operator O, defined in Eq. (|^), is the long distance approximation to the ly-exchange 
graph. The full amplitude will of course also involve a similar non-local matrix element but 
with the "singlet" operator O replaced by the octet operator Og- For now we concentrate on 
the singlet operator. None of the arguments given in this section depend on the particular 
choice of the operator. 

We will now argue that for heavy h and c quarks the non-local matrix element in Eq. (^ 
is well approximated by the matrix element of a sum of local operators. The approximation 
is valid provided Aqcd ^ "^c,b5 ^-c, the corrections are order Aqcd/^^c,?)- There are also 
corrections of order AqcD'^b.c/q'^- So our results are limited to the region were q^ scales hke 
m^j. The region were q^ does not scale like ?ti^^ is parametrically small, so the arguments 
we present are theoretically sound. However, there is the practical issue of determining 




u 



FIG. 4. Same as Fig. ^ but with the electromagnetic current couphng to the u-quark. 

a minimuni q^ for realistic calculations were our approximations can still be trusted. We 
return to this practical matter below, when we attempt to estimate the rate for this decay. 

The underlying decay is represented in the quark diagrams of Figs. 0-|^. In the heavy 
quark limit, Aqcd *^ ^^c,bi the heavy meson momentum is predominantly the heavy quark's. 
This suggests the following kinematics in the quark diagrams: for the momenta of the heavy 
quarks take mbV + ki, and tjIcv' + fee, for the momenta of the light quarks take ku and kd 
and then the photon's momentum is determined by conservation, q = rribV — irtcv' + Yl,ki. 
We can now exhibit our OPE by considering the quark Green functions in Figs. ^^. The 
convergence of the expansion for physical matrix elements rests on the intuitive fact that 
the residual momenta ki will be of order Aqcd (parametrically all we need is that these are 
independent of the large masses). This intuition is made explicit in Heavy Quark Effective 
Theory (HQET): there are no heavy masses in the HQET-Lagrangian so the only relevant 
dynamical scale is Aqcd- Thus our expansion of a non-local product will be in terms of local 
operators of the HQET. 

Calculating the Feynman diagram of Fig. with our choice of kinematics we have 

Here Qc = 2/3 is the charge of the c-quark and the tensor product corresponds to the two 
fermion bilinears. External legs are amputated. Using q = nihV — rricv' + J2ki and expanding 
in ki/uic^b we obtain, to leading order 
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FIG. 5. Same as Fig. |2| but with the electromagnetic current coupling to the d-quark. 

m^ — m^ 

This Green's function is that of a local operator in the HQET. Denoting by h^'' the anni- 
hilation operator for the heavy quark with four- velocity v, we define 

o = Jr, /^w h^:}r^ u. (7) 

Here V^^c are arbitrary Dirac matrices. With V^ ® Ff, set equal to the tensor product in (^, 

r. ® F, = Q.r'^^i^p. ® 7.P-, (8) 

mi — m^ 

the operator expansion is 

/ d^x e'^-" r[c7^c(x) 0(0)] = O + ■ ■ ■ . (9) 

The ellipses indicate terms of higher order in our expansion, and correspond to higher 
derivative operators suppressed by powers of wic,;,. There are also perturbative corrections 
to this expression. These show up as modifications to the operator defined by setting Fc(8>Fb 
equal to (^. 

The diagram of Fig. ^ can be analyzed in complete analogy. It leads to the operator O 
with the choice 

Te ® r, = -Q,7^P_ ® 7.P- ^'' + ^f 7^ (10) 



where Qb = —1/3 is the b quark charge. 

The analysis of Figs. ^ and ^ is similar, but there is an important distinction. With 
the electromagnetic current coupling to the light quarks, we get intermediate light quark 
propagators. The denominator in these propagators are parametrically large only if q"^ is 
parametrically large, i.e., if q^ ~ ^b- With this caveat, the OPE for Fig. HI gives 



Fe ® Ffe = -QuYP- ^1^ ® l.P- (11) 

r 



and for Fig. | the OPE gives 



Fe ® Ffe = Qdl^P- ® Y^luP- (12) 



III. SPIN SYMMETRY 



We have shown how to replace the time ordered product in Eq. (^) by a local operator. 
The replacement is valid provided the invariant mass of the lepton pair is large, i.e., scales 
as q^ ~ ^b- The operator O that replaces the time ordered product is defined by Eq. (^, 
with the tensor Fc (S> F^ defined as the sum of the contributions in Eqs. (H), (p!oD , (p!I|) 
and (0). We now show how to relate the matrix element of this operator to the operator 
with Fc ® Ffc = 7^P_ 'juP-- This operator is not only simpler, but one can estimate its 
matrix elements by a variety of means, as we explain below. 

Consider the matrix element of O as defined in Eq. (0) for arbitrary tensor product 
Fc (S> Ff, between heavy meson states. We will use heavy quark spin symmetry to determine 
the matrix elements of this operator between heavy meson states. Recall that the HQET 
lagrangian 

CuQET = h^'^tv ■ Dh^^^ + h^)%v' ■ Dh'^) (13) 

is symmetric under the group SU{2)b x SU{2)c of transformations acting on spin indices of 
the heavy quark fields: 



8 



"'V '-'b"'v ) 'V "-"c'V • 

At f ' = f the symmetry is enlarged to f/(4), which contains an SU{2) subgroup correspond- 
ing to a flavor symmetry. For now we will need only the spin symmetries. 

In order to make use of these symmetries, it is convenient to represent a spin multiplet 
consisting of a pseudoscalar P and a vector meson V^hj a 4 x 4 matrix 

^.= (^) [^.7^-^75]. (14) 

Then Sb G SU{2)b and Sc G SU{2)c act simply on the left, 

Hi''^ ^ SbHi'^ Hi^^ ^ SMi:}\ (15) 

while an arbitrary rotation R represented by the Dirac matrix 'D{R) acts simultaneously on 
both multiplets according to 

HiQ) -^ V{R)^H^^^V{R). (16) 

Consider now the matrix element {Hl,\0\HJJ'^). It must be linear in the tensors T^ ® Tf,, 
H^^') and H^j\ Acting with SU{2)b we see that Tj, -^ Tf,Sl and H^^^ -^ SbH^^\ so they enter 
the matrix element as the product VbH^\ A similar argument with SU{2)c gives then 

{H^\d\Hi'^) oc hI^^T, ® TbHi'^ (17) 

Finally, invariance under rotations implies that the remaining four indices must be con- 
tracted. There are two possible contractions, 

Tr(iTif)re)Tr(r5ij(^)) and TT{H^J^T,TbHi'^). (18) 

We now show that the second one is excluded by chiral symmetry. The lagrangian for a 
massless quark in QCD, 

C = ^iptp, (19) 

is invariant under the chiral symmetry 



^ ^ e*"'^^?/', (20) 

where a, the parameter of the transformation, is a real number. Under this symmetry the 
transformation rule for our tensors is 

TbH^''^ -^ e''''"'^rbHi''^e'""'' (21) 

and 

if^f^r, ^ e'"^«^ff^r,e-*"^^ (22) 



It is seen that the first contraction of indices in (|T8|) is invariant, but the second one is not. 
We have shown that heavy quark spin symmetry, rotations and light quark chiral sym- 
metry combine to give 

(i7(f)|d|iff ) = i/5(t.)Tr(^(f)r,)Tr(r,i7f)). (23) 

We have indicated that the invariant matrix element /3 is a function oi w = v ■ v'. In 
general, it is a function of v and v'. However, since it must be Lorenz invariant and since 
v'^ = v'"^ = 1, it is a function oi w = v ■ v' only. 

The octet operator in the HQET, 

Os = dT, T'^h^^ h^)v^ T'^u, (24) 

has the same spin and heavy flavor symmetry properties as its singlet counterpart. Therefore 
in complete analogy we can introduce a reduced matrix element /Jg: 

(ffif^lOsliff ) = ^/38(«;)Tr(^if)r,)Tr(r,i7W). (25) 

The authors of Ref. proposed a relation analogous to Eq. (^) for a AB = 2 transition. 
It was noted there that spin symmetry allowed more than one invariant and that, however, 
all invariants lead to the same symmetry relations. One may wonder if our use of chiral 
symmetry may help relate the different invariants there. We show that this is not the case. 
For the AB = 2 case the analogue of Eq. ( |T7|) is 
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(iff |OAs=2|i^f ) oc r-,^f ® nHi'\ (26) 

where Oab=2 = tFr^^ h^^'' JF^ /;,[,'') (note that we define h^^ to create a 6-antiquark). Again, 
invariance under rotations imphes that the remaining four indices must be contracted and, 
again, there are two possible contractions, 

Tr(r-,i?®)Tr(r5iff) and Tr(r^i?fr,/7W). (27) 

Chiral symmetry for the antiquark's meson tensor is just as for the quark's in Eq. (|2T|) , 

Tgi/f) -^ e-^"^T-bi?f)e^°^^ (28) 

Therefore both contractions in (^) are allowed by chiral symmetry. However, it is easy to 
see that for a class of operators of interest the two contractions are equivalent. If 



or 



for any arbitrary Dirac matrix T the two contractions are related by Fierz rearrangement. 
This class of operators includes the B — B mixing case studied in Ref. 0]. 

IV. QCD CORRECTIONS 

Consider the operator expansion in Eq. (|^). We have seen that at leading order the 
operator on the right hand side is given by Eqs. (0)^(il)- We now consider the leading-log 
corrections to this relation. In the large mass limit these are formally the largest, leading 
corrections to the operator expansion. A renormalization scale /x must be stipulated for the 
evaluation of matrix elements of the composite operators on both sides of Eq. (^. It is often 
convenient to evaluate the matrix elements at a low renormalization point /i = /iiow This 
choice makes the matrix elements in the HQET completely independent of the large masses 
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of the heavy quarks. If /iiow ^ ^^0,6 there are large corrections to Eq. (^ in the form of powers 
of as ln(mc,b/yUiow)- These powers of large logarithms can be summed using renormalization 
group techniques. The corrections to these "leading-logs" are of order 1/ ln(mc,b//iiow) or 
as. It is important therefore to keep /iiow small, but large enough that perturbation theory 
remains valid. When we estimate decay rates below, we use /iiow = 1-0 GeV. 

To study the dependence on the renormalization point /i we take a logarithmic derivative 
on both sides of Eq. (|^). Consider first the left side. Acting with fi{d/dfi) on the charm 
number current C'j^c gives zero, because the current is conserved. The action of fi{d/dfi) on 
the composite four-quark operator is a linear combination of itself and the octet operator. 
It is therefore convenient to consider instead the linear combination that appears in the 
effective Hamiltonian (^: 

d^x e''?-^ T[cYc{x){cO{0) + CsOsiO))] 

cO + CsOs + ■■■. (29) 

The coefficients c and Cg are such that the left hand side is //-independent. This is necessary 
for the physical amplitude to be independent of the arbitrary choice of renormalization point 
fi. Therefore our task is to determine the proper /i-dependence for c and cg so that the right 
hand side is also independent of /i. Therefore, if the operators satisfy 



_d_ 

dfi 






1 






(30) 



where 7 is a 2 x 2 matrix of anomalous dimensions, then the coefficients must satisfy 



/^ 



d 
dfi 



f A 



ycgy 



-7 



f~c\ 



ycsj 



(31) 



Here "T" denotes transpose of a matrix. 

The calculation of the anomalous dimension matrix is straightforward. In dimensional 
regularization with D = A — e dimensions, one needs 0] the residues of the e-poles of graphs 
with one insertion of the operators O and Og- The leading- log corrections arise from the 
leading, 0{as) terms in 7. These arise from the one-loop graphs in Fig. ^ 
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FIG. 6. One loop Feynman diagrams for the calculation of the anomalous dimension matrix. 
The solid diamond represents the local operators O or Og. 
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In principle the different tensor structures Tc ® Th defining O and Og can have different 
anomalous dimensions and even mix among themselves. However spin symmetry ensures 
that the anomalous dimension matrix is independent of the tensor structure Tc (S> Tfe. 

We find 



where 



/ 



7 



An 



8 -Awr{w) - 2 

4 17 14 



\ 



\^ -^wr{w) - 9 y - ^wr{w)) j 



(32) 



r[w] 



: ln(u7 + y/w"^ — 1). 



(33) 



The solution to the renormalization group equation (|3l|) is straightforward. In terms of the 
ratio of running coupling constants 

_ ( Ols{lj) 



and the functions 



MsM, 



1 41 - lAwr{w) 

12 ho 

3 1 + 2wr{w) 
^~'l h, ' 



(34) 



i, 



(35) 
(36) 



where the coefficient of the one loop term of the /^-function for QCD is 6o = H ^ f'^/; and 



Hf is the number of light flavors (rij = 3 in our case), we obtain 



( cif^) ^ 



V 



C8(/i) 



u 



c(/^o) ^ 



csif^o) 



(37) 



/ 



where 



U = z^ 



/ i^« + |^-? ±(,i-,-C 



27 



Z^ — Z 



l(A-y-li\ 8^€. 1^-? 



[z"- - z 



z^ + ^z- 



\ 



J 



(38) 



The question that remains is how to determine the coefficients c and Cg at some scale 
/Xq- But we have already determined these coefficients in Sec. |I|. Recall that the operator 
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O that replaces the time ordered product is defined by Eq. (^, with the tensor Tc ® r^ 
defined as the sum of the contributions in Eqs. (§), (|10D , (pJ]) and (|12]) with unit coefficient. 
The question can be rephrased as what is the scale /ig for which the calculation in Sec. |I| is 
valid. What we would like to do is to determine for what choice of /ig the loop corrections to 
relations like Eq. (^ will be free from large logs. The only relevant scales in the problem are 
the large masses rric^b, the invariant mass of the e^e~ pair, g^, which itself scales like m^^, 
the small masses and residual momenta and the renormalization point /ig. The corrections 
to the relations of Sec. || are guaranteed to be free from logs of the small masses or residual 
momenta. But there will be logs of ratios of large masses to the renormalization point, 
ln(mc,b//io). To avoid these one may choose /ig ~ rnc,b- For our computations below we 
will use /ig ~ 4.0 GeV. If the scales rric and nib are both large but very disparate one could 
review the above analysis by introducing a new renormalization group equation to re-sum 
the logs of mc/mb. The results of this section would still re-sum the logs of fi/rric. 
We thus have that c(/i) and C8(/i) are given by Eqs. (|37|) and (pSf), with 



2 1 

c(/^o) = c(/io) = 3(3;"^ - 2^^) (39) 

C8(/^o) = csil^o) = ^'^ + ^^ (40) 

where 

For illustration we have given the leading log expression for the coefficients c(/io) and 
C8(/ig), but in rate computations below we use the next to leading log results from 0. We 
do not have at present a full next to leading log result: still missing is a computation of the 
one loop corrections to the coefficients c and Cg at /i = /ig and of the anomalous dimensions 
matrix 7 of Eq. ( pOD at two loops. It is interesting to note that the coefficients c(/ig) are 
significantly enhanced at next to leading log order. For the case /ig = 4.0 GeV one has in 
next to leading order [Q c = 0.16, rather than the leading log result c = 0.07. We emphasize 
that this enhancement can be systematically accounted for. The large enhancement is not 
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a signal of perturbation theory breaking down but rather due to the accidental cancellation 
in the leading order. 

V. RATES: B^ -^ D'^*')^E+E- 

We are ready to compute decay rates. Defining 

/,Wm = ^D^*)o\jd^^ e^.- Ta^^(x)Kff(0))|5°), (42) 

the decay rate for 5° -^ D^*^^e~^e~ is given in terms of q^ and t = {pd +Pe+)^ = {pb —Pe-Y 
by 

2 



dV 






(43) 



dq'^dt 287r3M| 

where i'^ = u{pe-)l^'v{Pe+) is the leptons' electromagnetic current. A sum over final state 
lepton helicities, and polarizations in the D* case, is implicit. 

To compute h'^*'>^ we need to pull together the results of the previous sections. First 
the time ordered product is expanded in terms of local operators as in Eqs. (p|)-(p^. This 
involves replacing the coefficient functions c(/io) and c^^Hq) by c{j1q) and CsifJ^o) as seen 
in Eq. (|29|) . Then the matrix elements of the leading local operators O and Os between 
particular states can all be expressed in terms of the reduced matrix elements (3 and jSg 
defined in ( p3|) and (|25|) . Finally, to make all dependence on the heavy quark masses explicit, 
we run down the coefficients c and cg from the scale /i = /ig of order of rrib^c (which we take 
to be ^Jmjrii) to a scale /i = /^low of order of a few times Aqcd- 

Our computation gives 

K r — (2wmfc + mc)v^ — {rrib — Awm^^v'^ 



/i^ 



{mhV — rricV'Y 

3(mfef '^ + rricV^) 



ml — ml 



(44) 



and 
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K |-m;,(e'^ + 2v ■ ev^) — ■mc{3v ■ ev''^ + we'^) 



3 



{mi)V — rricV') 



vV'aV^ 



H "tAJl (45) 

(mfef — rricV'y 



ml — ml 

Here n = GF/V^V^iy*j^[cP + CsPs]- These expressions are our central results, demonstrating 
that the decay rates for B^ -^ D^*^^e^e'^ can be expressed in terms of the matrix elements /? 
and /?§. Below we make an educated guess of these matrix elements, but for reliable results 
they should be determined from first principles, say, by Monte Carlo simulations of lattice 
QCD. 

In the computation of the rate the amplitude depends on heavy quark masses m^c and m,/,, 
while the phase space involves physical meson masses Mb and Md or M^.. Although it is 
straightforward to retain the dependence on all four masses in our expressions for the decay 
rates, we have chosen to express the results in terms of physical meson masses, with the 
substitutions m?, = Mb and m,c = Md or m,c = Md*. We are not justified in distinguishing 
between quark and meson masses since the distinction enters at higher order in the l/m,c,b 
expansion. 

It is now a trivial exercise to compute the differential decay rate. Integrating the rate in 
Eq. (E3|) over the variable t we obtain 



dV a'G 



2r^2 



V,bVud\\c(3 + c^f3syT{q). (46) 



t/g2 2887rM| ' 

Here J^{q) is a dimensionless function of g = \jq^lmt^ and m = Mj^(*)/Mb- For B^ 

D*^e~^e~ it is given by 

4 VI - 2g2 - 2m^ + q^ - 2m^q^ + m^ 
3 q^m{l — m'^y 

{5m^ + 19m^g2 + 30m^ - 20m^ - Ufm^ 
-20m^ + 12 fm^ + f + f- 2q^ + 2m^q^ 
—Qm^cy + 5m^° — Q(fm} + 5m^q^), 



(47) 
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while for B^ -^ D^e+e' 

4(2m2 + 1)2(1 - 2f - 2m2 + q^ - 2m'^f + m^)i 

3g4m(l-m2)2 • ( ^^ 

In these we have neglected the electron mass. 

In order to obtain a numerical estimate of the branching fraction we need to calculate the 
hadronic matrix elements (3 and /Sg- While these could be studied in Monte Carlo simulations 
of QCD on the lattice, at the moment we have no reliable information on their magnitude. 
These matrix elements are similar to the matrix element of the AB = 2 operator ioi B — B 
mixing. Lattice QCD |^ indicates that the vacuum saturation approximation works very 
well for B — B mixing. Therefore we take vacuum saturation as an educated gues^ for f3 
and Ps- Taking Fc (S> T;, = •y^'y^ <S) 7*^75 the right hand side of Eq. (^) is v^v''^(3. On the left 
hand side vacuum saturation gives {z~'^' f bV'b I ^ ^b){.z~'^' f dp'^q / ^/ Mjj) . Here z is defined 
in Eq. (P^ and aj = 2/5o [0 is the well known anomalous scaling power for the heavy-light 
current in HQET.^ Thus we obtain 



(3(w) = z-^^'fBfD^MBMo (49) 

/38H=0 (50) 

The second equation is true not just in vacuum saturation but also in the approximation 
that we can insert a complete set of states between the currents defining Og. This is not 
an exact statement because the composite operator O^ does not equal the product of two 
currents. But the distinction arises from their different short distance behavior. So we 



^The matrix elements in B^ — > D^*''^e^e can be related by symmetry to the matrix element for 
B — B mixing, if the matrix element of the octet is negligible; see Ref. |^] 

^The two factors of z~°'' really correspond to distinct running, between tti;, and //jow for the first 
factor, and between rric and ^uiow for the second. The distinction is of higher order than we have 
retained, if we assume that the heavy scales rrih and rUc are not too disparate, that is, that Us does 
not run much between these scales. 



expect the deviation of Ps from zero to be of order of the QCD couphng at short distances 
a(/io) times the unsuppressed (3. 

Using these matrix elements we integrate the differential rate in Eq. p6| ) over the range 
1.0 GeV < q^ < g^ax ^o obtain a partial decay rate. We have chosen g^j^ = 1.0 GeV as a 
lower limit since our OPE requires that g^ scale like m'^fj. The corrections to the leading 
terms in Eqs. (0) and (0) are of the form of an expansion in rric^bk/q'^, where k is any of the 
residual momenta and in our matrix elements is of order Aqcd- Parametrically, if g^ ~ m^ j, 
then rric^bk/q'^ ~ AQco/^^c.b <^ 1. In addition, the region over which g^ < Aqcd"^c,6 where 
the expansion breaks down, is parametrically small. However, physical heavy masses are 
not very large, and the scale mftAqcD is just slightly smaller than m^. In order to have some 
non-trivial phase space we have taken g^ > mfeAqcD ~ 1.0 GeV. The price we pay is that 
for the lower values of g^ our expansion converges slowly, nic^bk/q^ ^ 1. 

We find 

Br(S° ^ D*°e+e-)|,2>i cev = 1-4 x 10"^ (51) 

Br(5° -^ L)°e+e-)|,2>i cev = 2.6 x 10"^ (52) 



where we have used iKfeKidl = 0.04, /£> = JeJ Mb/Md and fs = 170 MeV. It is important 
to observe that the portion of phase space g^ > 1.0 GeV is expected to give a small fraction 
of the total rate since the pole at g^ = dramatically amplifies the rate for small g^. The 
rates for 5^ -^ D*^e^e^ and B^ — > D^e^e^ can be obtained to good approximation by 
replacing |V^6\4d| by iKbKisI, reducing the rates by (0.22)^ ^ 0.05. 

The next generation of B-physics experiments at high energy and luminosity hadron 
colliders, like LHC-B and BTeV, will produce well in excess of 10^^ S-mesons per year. Our 
calculation includes only large invariant mass lepton pairs so detection and triggering on 
the lepton pair should be straightforward. Dedicated studies must be done to determine 
feasibility of detection and measurement of spectra of these decays. 
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FIG. 7. Feynman diagram representing a contribution to the Green function. The filled square 
represents the four quark operator O and the cross represents the electromagnetic current jj^^, cf. 
Eq. (|58[), which here couples to the c-quark. 

VI. DECAYS TO QUARKONIUM 
A. Operator Expansion and NRQCD 

The decays Bs -^ rice~^e~ and Bs -^ J/ipe'^e^ (and obvious extensions to excited char- 
monium) can be studied in a similar way. The notable difference in the operator expansion 
here is that the residual momenta k of the heavy quarks in the quarkonium bound state 
do scale with the large heavy mass k ~ asfric, as opposed to the residual momenta of the 
quarks in the heavy B ot D mesons, k ~ Aqcd- The residual momentum for the case of 
quarkonia is small for a different reason: k = rricU and k^ = ^rricU^ are small because the 
velocity u of the bound quarks is small for heavy quarks, u ~ as{mc). The parameter of 
the expansion is therefore mc^bk/mli^ ~ asirric). 

Our best hope in making the nature of the expansion explicit is to use NRQCD [Q, 
the effective theory of non-relativistic quarks in QCD. As opposed to HQET, where all the 
heavy mass dependence has disappeared, the lagrangian of NRQCD still depends on the 
heavy mass: 

jC^rqcd = -^HtDt - ^)* (53) 

Here \1/ denotes a two component spinor field for the c-quark. A separate spinor field 
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must be included to describe the antiquark. We have written the lagrangian in the rest- 
frame of charmonium, but it is straightforward to boost into a moving frame. One rehes 
on the dynamics to generate the small parameter of the expansion.^ For example, the 
two terms in £nrqcd are of comparable magnitude if, as expected, Dj ~ A;° ~ nical and 
|D| ~ |k| ~ rricas. 

The operator expansion is in terms of operators with an HQET quark, a light quark and 
a pair of NRQCD quark- antiquark. So instead of Eq. (|^) we have 

O ^ dT, /if vl/tr^ vl/,, (54) 

where \1/J and \l/c create a charm quark and a charm antiquark, respectively. We elect to use 
four component spinors throughout; the reduction to two components results from algebraic 
constraints that must be imposed, just as in HQET: 

The calculation proceeds much as before. The effective Hamiltonian for the weak tran- 
sition is 

Kff = -^ K.K; ic{^^/Mw)0 + Cs{fi/Mw)Os) , (55) 

where 

= sYP-h c-f^P-C (56) 

and 

Os = sYP-T% c-f^P-T^c. (57) 



^Attempts to make the expansion in u Q or, alternatively, in 1/c [ p!o[ | explicit yield theories where 
the gluon self-couplings must be perturbative. The scale of QCD must then be negligible compared 
with the Bohr radius of quarkonium, Aqcd ^ rncas{mc). In our case non-perturbative gluons play 
a crucial role in binding the heavy-light meson B. 
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FIG. 8. Same as Fig. but with the electromagnetic current coupling to the c-antiquark. 




FIG. 9. Same as Fig. ^ but with the electromagnetic current coupling to the 6-quark. 
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The operator expansion of the hadronic matrix element takes the form 

d'x e^^- T[ji:Jx){cO{0) + CsOsiO))] = 

cO + CsOs + ■ ■ ■ , (58) 



where O is defined in (B^ and the octet operator Os is defined analogously, 

Os = dT^ T'^h^^'^ ^Jr, T-^^g- (59) 

The first task is to determine the tensor F^ ® T^. To this order we consider Green functions 
of the time ordered product in Eq. (|58D with four external quarks. The in-going momenta of 
the h- and s-quarks are mhV + fcb and fc^, respectively. The outgoing momenta of the charm 
pair are ttIcv' + k^ and nicv' + kc- As explained above, we expect kf, r^ kg r^ Aqcd while 
kc ^ kc ^ ?7ic«s('^c)- The leading term in the momentum of the electromagnetic current is 
q = rUhV — 2771 cv' . For the purpose of determining the expansion coefficients at tree level we 
may set c = 1 and Cs = and, choosing a renormalization point /ig of the order of the large 
masses 77ic,fe, we can set c = 1 and cs = 0. There are four graphs contributing to the tensor 
Te ® Tfo. Fig. gives 



^7 7lbi) -^cif - 1 ) 
77ll — 277lb77lcW 



Tc ® r, = Q.r ''''^~T::s^:^ YP- ® i.p-, m 



and Fig. |] gives 



r.or> = ftrP- "'"f^:"'''''^" v«7.P-. (61) 

771^ — l77lij77l(.W 

Note that the denominator, which dictates the convergence of the expansion, scales with 
m^j. It vanishes at Wq = 7nf,/277ic. However, this is never in the physical region: Wmax = 
(m^ + 4:77il)/4:77ib77ic = wq — {77ib/4:7nc — 77ic/77ib), but 771^ > 2771c for the decay to be allowed. 
The diagrams in Figs. ^ and |1^ are just as in Figs. ^ and ^, with the replacement 
q = 77ibV — 77icv' —^q = rribV — 2mcf '. For the first we have 

Te ® r, = -QbYP- ® luP- '^'.^^Ti ^"^ (62) 

m^ — 4m^ 
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FIG. 10. Same as Fig. |^ but with the electromagnetic current coupling to the s-quark. 
and for the second 

Tc ® r, = g,7^p_ ® r^iuP- (63) 

r 

Again we see that the expansion remains vahd as long as q^ scales with the heavy masses 
(squared), and this limitation arises solely from the coupling of the photon to the light 
quark. 

B. Spin Symmetry 

The NRQCD lagrangian contains separate fields for the charm quark and antiquark. The 
quark lagrangian, Eq. (|53D is symmetric under spin-5'f/(2) transformations. The antiquark 
lagrangian is similarly invariant under a separate spin-5'f/(2). This case has a larger spin 
symmetry than the case of decays to D-mesons. One can therefore write a trace formula 
analogous to Eq. (p3D without using chiral symmetry of the light quarks. 

We can represent the charmonium spin multiplet (?7c, J/ip) by the 4x4 matrix 

i,»)=(i±i:)|Vvr-..7d(^). (64) 

The action of spin-S'[/(2) x SU{2) on this is then 

H^f -. SMifsl (65) 
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Consider the matrix element {Hi? \0\H^^''). It must be linear in the tensors FcC^rb, if^''^ 
and Hyf . As before, acting with SU{2)b we see that Ff, -^ ^bSl and H^^^ -^ SbH,ll'\ so they 
enter the matrix element as Vi,H^\ Now, acting with the spin symmetries of NRQCD, we 
have Eq. (^) and T^ -^ ScTcSl, so that they must enter the matrix element as Tr(^^fTc). 
Finally, rotations demand that we sum over the two remaining indices, 

{HiflOlHi"^) = ^/3Tr(^if)F,)Tr(F,iff ). (66) 

Similarly, for the octet operator we find 

{Hif\ds\Hi'^) = ^Ps{w)TT{HifT,)TT{TbHi'^). (67) 

We have used the same symbols here for operators and reduced matrix elements as in Sees. H 



and III, but they should be understood as distinct. 



C. QCD Corrections 

Consider the operator expansion (pH]). Just as in Sec. |T^ we argue that matching between 
left and right sides is most conveniently performed when the renormalization point /ig is 
chosen to be of the order of the scale of the heavy quarks. For simplicity we assume that 
rric and rrih are not too different, but very big, so that we do not have to worry about large 
logs of the ratio mc/mb. Then one may take, say, /ig ~ ^A'V"^- The point is that the 
coefficients on the left hand side of ( ^8|) explicitly depend on Miy//io and the operators 
implicitly depend on rricfi/ ^Iq. If we choose to do the matching at a scale /Xg ^^at differs 
much from rrtc^b then there are implicit large corrections. Note that the right hand side 
of (|58|) can only introduce logs of low scales over /ig, but the same infrared logs are found 
on the left side of the equation. 

Once the coefficients c and cs in (^8]) have been determined at /iq we must ask at what 
scale /i we should evaluate the matrix elements and how to get there. The situation is 
more complicated than in the case of B^ -^ D^e^e^ of Sec. |V| because now the matrix 
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element in the combined HQET/NRQCD effective theory has several scales. In NRQCD 
the relevant distance scale is the inverse Bohr radius mcO-si'mc) and the relevant temporal 
scale is the Rydberg mcolimc). In HQET the dynamical scale is Aqcd- Of course Aqcd also 
plays a dynamical role in NRQCD, but it is usually taken to be irrelevant since one assumes 
Aqcd ^ 'mcoilirric) -C rricasimc). So we are faced with a multiple scales problem. Setting ix 
equal to any one of these scales leaves us with large logs of the ratios of /i to the other two. 
It is not known how to use the renormalization group equation to re-sum these logs. 

Suppose that we set /x ~ mca{mc) or yU ~ mca'^{'mc). If we then use the renormalization 
group to sum powers of a;(mc) ln(r7ic//i) we will be summing powers of a(rric) In a (mc). 
Notice that these logs vanish as m^ -^ oo, since a(mc) ~ l/ln(mc/AQCD)- Contrast this 
with the case /i ~ Aqcd (or, generally, setting /x equal to any fixed scale as nic — > oo). Then 
a(mc) ln(mc//i) ~ 1 as rric — * oo. As a matter of principle, in the large mass limit it is 
these latter logs that must be summed (they are parametrically of leading order in the large 
mass expansion). Therefore we re-sum the leading logs with a fixed low scale /x = /xiow and 
choose, as before, /iiow = 1-0 GeV in our numerical computations. 

In order to use dimensional regularization and keep track of different orders in the non- 
relativistic expansion we adopt the 1/c counting advocated in Ref. |[TD|. However, we use 
a covariant gauge for our calculations. This is convenient because the Feynman diagrams 
involve light and HQET quarks in addition to the NRQCD quarks. In leading order in the 
1/c expansion the quark lagrangian in (^3]) is replaced by 

>Cnrqcd -^ ^^(^A - 7^)^ (68) 

znic 

The only interactions are due to temporal gluon exchange. Since we work in covariant gauge, 
this is not a pure Coulomb potential gluon. It is easy to see that no diagram involving an 
NRQCD quark gives a divergent contribution. The self-energy diagrams for the NRQCD 
quarks have an infinite piece, which however is independent of the momentum and therefore 
gives no contribution to wavefunction renormalization. Therefore the four quark operators 
scale as the heavy-light currents. That is 
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f.n\ 



4 



v°i/ 



(69) 



is the anomalous dimension matrix in the renormahzation group equation for the operators, 

f n\ f n\ 

(70) 



^^ 



d 

dfi 



O 



1 



o 

yd, J 



Then the coefficients must satisfy 



/i 






f,\ 



-7 



(71) 



ycgy 



where, as above, "T" denotes transpose of a matrix. 
The solution is trivial, 

c{fi) = z'^'cifio) (72) 

C8(/i) = z-^'^'csifio), (73) 

where z is defined in Eq. ( p4D and aj = 2/bo is the well known anomalous scaling power for 
the heavy-light current in HQET ||^. 

Contributions from higher orders in the 1/c expansion produce mixing with higher di- 
mension operators and are therefore excluded to the order we are working. This is easy to 
see. To compensate for the powers of 1/c one must have additional velocities in the oper- 
ators. But these come from powers of d/rric. The leading correction to the lagrangian is 
of order 1/c'^/^. Since two insertions are needed this gives a graph of order 1/c^. Since one 
power of c is needed to form the QCD fine-structure constant, a^ = g^/inc, the divergent 
part of the graph involves p'^/mlc^. It is straightforward to verify this by direct calculation. 



D. Rates 



Defining 



/^Wm = (vl;| / d'a: e^^^ T{jUx)KAmBs), 



(74) 
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where ^ = t]c, J/ip, the decay rate for Bg — > \l/e"'"e is given in terms of g^ and t 
(p* +Pe+f = {pb -Pe-f by 

9 2 



dr 



1 



^,,/iW'^ 



2^M 



(75) 



rfg^cft 2%3M3 q 

where C^ = u{pe-)'y^v{pe+) is the leptons' electromagnetic current. A sum over final state 
lepton helicities, and polarizations in the \E' = J/ip case, is implicit. 
We obtain 



}^{Vc)fJ. 



K vmiyv'^ — 2(wmb — mc)v^ 
3 L {miyV — 2mcV'Y 



+ 



mt 



Ami 



(76) 



and 



h{J/i')t^ 



K r2mhV ■ ev^ — {nib — 2mcw)t^ — 2mcV ■ ev'^ 



l\2 



+ 



{nihV — 2mcV 



[mfyV — 2'mcV'y ' ml — 2ini,inc.w 
m^e^ + 2mc{v ■ ev'^ — we^) + 2imce^°'^'^eaVpv'^ 
ml — AmA 



(77) 



Here k = GF/V^V^jy*g[cP + Cg/^s]- These expressions are our central results for decays to 
charmonium, demonstrating that the decay rates for Bg -^ T]ce^e~ and Bg — ^ J/ipe~^e~ can 
be expressed in terms of the local operator matrix elements (3 and jSs- 

We now compute the differential decay rate. We integrate the rate in Eq. (^) over the 
variable t and obtain, for both Bg — > r]ce~^e~ and Bg -^ J/ipe'^e', 






"'^^ 'V,iycs\\c(3 + c,(3,fHq)- 



2887rM| 



(78) 



Here J^{q) is a dimensionless function of g = Jq^ Iml and rh = Mj/^/Mb- For Bg 

J/ipe~^e~ it is given by 



J^ 



A^/l - 2q^ - 2m^ + q^ - 2m^q^ + m^ 
3q^rh'^{l — m'^y{l + g^ — m^)^ 
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-lllm^g^ + 234m^g^ + 104m^g^ + 92m^g^ 

-124m^g^ + m^^ - 72m^g^ + 55m^°g^ - 12m^2^2 
+23m^g^ - 78m^f + Am^q^^ - 6m^q^^ - 48m'^f] 



(79) 



while for Bg —>■ rj^e^e 



4{l-2q^ -2m^ + q^-2m^f+m^)2 
~ 3g4m2(l - m2)2 " ^^°^ 

For a numerical estimate we need to calculate the matrix elements /3 and /^g. Again we 
use vacuum saturation. However, now this approximation is supported by NRQCD. It is 
argued in Ref. [^ that soft gluon exchange with the quarkonium is suppressed by powers 
of the relative velocity u = aslrric), and that the matrix element of the octet operator is 
similarly suppressed. Therefore we take 



f3{w) = z-'^'fBU^MBM^^ (81) 

P8{w)=0. (82) 

Note that because vacuum saturation here is valid at least as a leading approximation in 
a velocity expansion, the combination of coefficients in ([721)- (ff^) and matrix elements in 
(pl|)-(|8^) is automatically independent of the renormalization point fi. Spin symmetry gives 
fvc = fj/ip- We use the measured value from the leptonic width in the tree level rate equation. 



. 2 ■/'j/V' 



r{J/^ -. e+e-) = Ana'i^, (83) 



and obtain fj/^ = 0.16 GeV. 

Integrating over g^ > 1.0 GeV we have partial branching fractions 

Br(S, ^ J/^e+e-) |,2>i cev = 2.2 x lO^^^ (84) 

Br(5, ^ r^.e+e-) |,2>i cev = 3.4 x 10^" (85) 

29 



where we have used iKbKsl = 0.04, and /b = 170 MeV. Again, we remind the reader that 
the portion of phase space q^ > 1.0 GeV is a small fraction of the total rate since the 
pole at g^ = dramatically amplifies the rate for small g^. The rates for B^ -^ J /ipe^e' 
and -B° -^ rjcC^e^ can be obtained to good approximation by replacing |V^f,Vcs| by |V^f,V^d|, 
reducing the rates by (0.22)^ ^ 0.05. The rate (^^ may seem too small to be detectable 
even in the next generation of hadronic colliders. However it must be kept in mind that the 
signature involves four leptons with large invariant masses (one being the J/ip)- 

VII. CONCLUSIONS 

We have successfully shown how to implement the OPE advertised in Ref. ^ to the 
processes Bd,s -^ J/ipe'^e', Bd,s — ^ Vc^^^'y Bd^s -^ D*^e^e~ and Bd,s -^ D^e^e" . By 
the use of the OPE the long distance (first order weak and first order electromagnetic) 
interaction is replaced by a sum of local operators. The application of the OPE is restricted 
to a limited kinematic region. 

In the processes Bd^s -^ J l'^e^t~ and Bd^s — ^ ^ce'''e~ our method leads naturally to an 
NRQCD expansion for the J/?/; and ric- This illustrates that the methods of Ref. [|^ are 
applicable to a wider class of processes. 

Furthermore we found that the number of independent matrix elements of the local 
operators is severely restricted due to a combined use of heavy-spin, rotational and chiral 
symmetry. The independent matrix elements could be determined, say, in lattice simula- 
tions. Our paper shows that the processes considered can be studied in a systematic fashion 
independent of any model assumptions in the kinematic regime of g^ scaling like m^^. 

Using a crude estimation of the matrix elements, we found the rates of all the processes 
considered to be small. We expect some of them, in particular Bg —* D*^e^e^ , should be 
accessible at planned experiments at hadron colliders, like BTeV or LHC-B. 
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